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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM 1272 


CRITICAL VELOCITIES OF ULTRACENTRIFUGES* 
By V. I. Sokolov 


In chemical industry, it is often necessary to deal with the 
phenomenon of instability of operation of ultracentrifuges that lead, 
in a number of cases, to breakdown, accidents, and early wear of 
the' instrument. In these cases, understanding of the phenomenon that 
may help to remedy the condition is not always found. 

The characteristics of the passage through the critical velocities 
of ultracentrifuges are discussed and the corresponding practical con- 
clusions are presented in this paper. 

The Euler differential equations of rotation of a solid body are 
applied. 

The following constructions are made: the -fixed axes £, q, and 

£ with the origin at the point 0 (fig. l)j the axis of symmetry of the 
centrifuge rotor z', the nodal line I perpendicular to the plane of 
#, z; and finally, the axis k perpendicular to the axes z, I. 

By applying the Euler equations of the rotation of a solid body and 
according to the notation of Nikolai (reference 1) 


d0 l 

dt~ + q l a z ” r l °k = L I 


da k I 

W + r l °1 - p l a z “ L k > 


d0 z 

+ pi a k - qi a i = L z 

J 


(i) 


The rotation of the centrifuge rotor may be considered as a regular 
precession of the symmetrical gyroscope. The angular velocity cd of 
the rotor is therefore the angular velocity of this precession and the 
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axis $ is the axis of precession. By laying off the angular velocity 
of the drum, or more accurately, the trihedron formed hy the axes I, 
k, and z that are assumed constant on the axis and by projecting 
it on the axes I, k, and z, the following equations are obtained: 


Pq = 3 
a a) sin (3 




( 2 ) 


r^ = | o> + a 


cos p 




The velocity a is the result of the action of the Coriolis forces 
and, by comparison with cd, is a small magnitude. For simplification of 
the solution, it is assumed that 


r 


i 


ca cos 


3 


(2a) 


In other words, the assumption, often used in the theory of 
vibrations, is made that the principal axis of the rotor z always 
remains in the plane including the vertical (the initial position of 
the axis of the shaft) and rotates with the same constant angular velo- 
city oo with which the drum rotates. 

The center of inertia of the rotor is assumed to lie on its axes 
of symmetry 0q and 0g. The axes I, k, and z will then be the 
principal axes of inertia of the rotor at the point 0-^ because, from 

the symmetry in regard to the axis z, all straight lines drawn perpen- 
dicular to this axis through the point 0q are principal axes of 
inertia. From this assumption, it follows that the moments CJj, CJ-^, O z 
are 


0 1 - Ap 

a k = Aq 

a = Cr 
z 

where A and C are the equatorial (at point 0q) and polar moments 
of inertia. 
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From equations (2) 

0 ± = Ap 
0^ = Am sin (3 
O = Cm cos p 

z y 



By linearizing, the following equations are obtained: 


d 1 = A(3 
q., = mp 
O k = Amp j 



(4) 


In equations (l), L k = L z = 0 but the moment of the external 

forces applied to the rotor relative to the axis I is determined by 
the expression 


L 


1 



sin (3 + Gh n 


sin p 


(5) 


where 8-^ is the force applied at the end of the rotor spindle that 

produces unit rotation in the case where the spindle is immovably 
fixed, and G is the weight of the rotor; h and h-, are shown in 

figure 1. 


' By neglecting the gravitational force and by replacing sin P 

■by 3 



( 6 ) 


By substituting in the Euler equation the values of the magnitudes 
that are used, the follo\d.ng equation is obtained: 



P = - 5 


'll 


h 2 p 


(7) 


IL 



4 


NACA m 1272 


Because the nonuniform distribution of the mass of the rotor 
results in a deviation of the principal axis of inertia from the verti- 
cal by the angle 3 q, equation (7) must be written as 



A P + 3 0 



where 


( 8 ) 


3 = 3 0 



CD 


O) 


Hence the critical velocity is 


co, 


'cr 



( 10 ) 


In order to complete the solution and obtain a computation formula, 
S -| -| is determined.. 


In the corresponding constructions of the usual type of centrifuges, 
the spindle is assumed to be of conical shape. In determining the 

spindle may be considered as a beam with the large cross section built 
in. The fixing of the beam corresponds to the joining to the spindle 
with the centrifuge rotor (fig. 2). 


On the basis of the Mohr formula, the displacement of the end of 
the beam A-^ may be determined by the equation 


pcw-Z 


A 11 = 


EI 


dx 


(ID 


a 


where M is the bending moment, E is the modulus of elasticity, and 
I is the moment of inertia of the beam cross sectionj a and l are 
shown in figure 2. 


The moment of inertia at any section of the spindle is 



4 


(12) 
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where r is the radius of the spindle cross section. As follows 
from figure 2 , r = kx and therefore 

I = gkjxj; ( 13) 

4 


Conversely, the tending moment due to unit force applied at the 
distance a from the origin of the coordinates (at the tip of the 
team) is 


M = x - a 


The expressions for I and M are substituted in equation (11) 
to give 

Pb 

|x - a | 2 


11 itk 4 E 


dx 


(11a) 


By integrating and substituting the limits 


A- 


11 ' „k 4 E 


1 a 1 sr 

3a ^2 b ^3 


(14) 


By denoting the largest diameter of the spindle by R and the 
smallest diameter by r^, k, a, and b are determined. In accordance 
with figure 2 


k = 



(15) 


a 



0 


r Q* ^ 
R - r 0 = 



(16) 


b = 7 + 


l 


(17) 
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Finally 



4Z V 


jiE R 



+ 



1 

R 



( 18 ) 


The obtained result is substituted in the formula for on.,, to 

give 


O) 


cr ~i 


h 2 jtE R - r 0 | 3 

1 A - C| 

r 0 1 1 r o 2 

4Z 3 

r 2 3r 0 R 3R 3 



(19) 


Because ultracentrifuges operate at velocities higher than the 
critical, it is impossible to ignore the passing of the rotor through 
the critical point. 

Factors that affect the passage of the centrifuge rotor through 
the critical velocity are to be considered. 

Two methods exist by which the passage of the ultracentrifuge 
rotors through the critical velocities may be affected. The first 
method provides a rapid start that excludes an inadmissible increase 
in energy of the forced vibrations and may be applied to small centri- 
fuges of various types. The second method restricts the amount of 
bending of the spindle when passing through the critical point by means 
of rings formed by suitable plates. The second method is the one 
applied most often and will be subsequently considered with the aid 
of a method proposed by Kapitsa (reference 2). 

As has been shown by Kapitsa, the critical point separates two 
qualitatively different types of rotation. The transition from one 
motion to the other for the ultracentrifuge is discussed. 

In the first period of the motion, the gyroscopic couple balanced 
by the elastic force of the spindle is determined by the left side of 
equation (8). 

For a certain value of the angular velocity , the head of the 
rotor touches the restricting surface. The gap between the head and 
the ring is denoted by e. The angle of inclination of the rotor axis 
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corresponding to the contact of the restricting ring will he equal to 
p^ = e/h 3 , where hg is the distance from the lower support to the 
ring. 

The value of (3^ is substituted in equation (8) and is solved for 
by replacing -h^S j j /C-A with the value • Thus 


o>l = m cr 


e 

h 3 P 0 + e 


( 20 ) 


When the velocity of the rotor reaches the value the deflection 

of the shaft ceases to increase because the restricting ring prevents 
it. 


After contact and with increase in angular velocity, the head of 
the rotor presses against the restricting ring through the centrifugal 
force. At first the head slides along the surface. When the friction 
force between the head and the ring reaches a certain value, however, 
the rotor starts to roll along the ring in the opposite direction and 
the spindle immediately straightens out. 

The value of eng corresponding to the straightening out of the 
spindle is determined. 

The contact of the rotor with the ring changes the value of the 
moment of the external forces determined by equation (l). The moment 
of the reaction force of the ring Q.h -3 is added to the moment of the 
elastic force of the bending spindle. 

The angle between the plane of action of the gyroscopic moment and 
the bending plane including the point of contact is denoted by y; the 
moment acting in this case is determined by the equation 


or 


C-A 


0]_ + P 0 cos y 


* By denoting as Q the reaction force of the ring, the following 
equation may be written: 


or 


C-A 


P-, + p n cos r 


= - 5 h c p - Qh 3 


( 21 ) 
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The sliding of the head of the rotor along the ring is accompanied 
by the appearance of the tangential force pQ, where p is the coef- 
ficient of friction. The . condition for the existence of motion of the 
first kind must be 


Pq q)^ A - C sin y 

h^ 


>Qh 


( 22 ) 


From equations (21) and (22), the following equation is obtained 
by the Kapitsa method (reference 2): 


— sin x - cos x p 


Pi a^2 . 


ao. 


cr 




(23) 


^2 Pl 



(23a) 


where co, is the angular velocity corresponding to the straight- 
ening of the spindle. The condition of the passage through the 
critical velocity is obtained in the form 

(24) 

By denoting the gap between the head of the rotor and the ring 

as e 


_Z_> 

h 3P 0 



(25) 


The preceding inequality shows that the ultracentrifuge rotor cannot 
pass through the critical velocity if the space between the head and 
the ring is small, if the coefficient of friction between them is small, 
or finally, if the unbalance of the rotor is large. 
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The increase in the gap e, however, is limited by the strength 
of the spindle. 

The maximum permissible value of e is determined. The displace- 
ment of the end of the rotor spindle under the action of a unit force 
is The stress arising at the place Where the spindle is stopped, 

in the case where the force x is applied to the end of the spindle, 
is 

a = xl/w 


where W is the moment of resistance of the section and the deflection 
is equal to Anx. By giving a value of the deflection equal to e 

corresponding to the occurrence in the critical section of a stress 
equal to the permissible stress a z 


e = A-^x = 


A i;l ct z W 


In accordance with equation (18) 


,2 4 

7 R u. 


e = 


E R - r 


0 + _ 1 _ _ 1 
2 3r q R 


r 0 

3R 2 ' 


or 


7 2 a 


e = 


E R - r 0 


R 2 + R 


3r, 


3 r 0 2 ^ I 

- R - 4- 


(26) 


(27) 


(23) 


The pressure at the upper bearing and at the instant of the passage 
of the rotor through the critical velocity is determined. 


In equation (26), the pressure on the bearing is evidently equal 
to the magnitude x or 


x 


n z W 

l 


(29) 


After passing through bhe critical point, the spindle remains 
bent by an amount practically equal to PqR . The force due to the 

bending transmitted on the bearing is equal to 


Xi = 


>11 


(30) 
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or by substituting the value from equation (10) 

_ o > cr 2 |a - c| p 0 


(31) 


CONCLUSIONS 

Particular attention is to be paid to the problem of the passage 
of an utlracentrifuge rotor through the critical velocity. The dif- 
ficulties encountered in the practical use of ultracentrifuges are 
connected with this passage. 

The dynamical balance of the rotors has often been neglected, 
whereas their balance increases in connection with the too frequent 
removal of the rotors from the body of the centrifuge for cleaning. 

A result of the increase in the unbalance may be the loss in the 
ability to pass through the critical velocity and the rapid wear of 
the bearings (the force acting on the bearing, equation (31)). 

The preceding discussion points out the factors that affect the 
passage of rotors through the critical point, namely, the gap between 
the head and the ring, the coefficient of friction, and the unbalance. 
These factors cannot be excluded in the operation and design of 
ultracentrifuges . 


Translated by S. Reiss, 

National Advisory Committee 
for Aeronautics. 
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